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(57) Verfahren zur digitalen Interpolation von Signa- 
len durch Multiplikation verzOgerter, mit einem ersten 
Taktsignal (f1) verkoppelter Eingangswerte (x(k)) mit 
entsprechenden Gewichtungsfaktoren (Ci. C|). die von 
einem Zeitdifferenzwert (At; dt) abhangig sind. der durch 
den Interpolationszeitpunkt (tP) und das Zeitraster des 
ersten Taktsignals (f1) bestimmt wird. wobei die Gewich- 
tungsfaktoren (Ci; C t ) durch eine Impulsantwort h(t) im 
Zeitbereich bestimmt werden. Die zugehdrige Ubertra- 
gungsfunktion H(f) weist im Frequenzbereich ein Signal- 
dampfungsverhaiten auf. das bezuglich der 



Sperrbereiche im weserrtlichen auf die bei den Fre- 
quenzvielfachen des ersten Taktsignals (f1) liegenden 
Starbereiche (1S. 2S. ...pS. ...) beschrankt ist Jedem 
dieser StOrbereiche werden mindestens zwei nebenein- 
anderliegende Nullstetlen zugeordnet. Oder beim Vor- 
handensein von Nullstellen doppelter Ordnung wird 
mindestens einem der StOrbereiche (pS) und den zuge- 
hOrigen periodischen StSrbereichen (3pS. 5pS, 7pS, 
9pS. ...) mindenstens eine weitere Nullstelle der Uber- 
tragungsfunktion H(f) zugeordnet. 
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EP 0 696 848 A1 

Beschreibung 



Die Erfindung betrrfft ein Verfahren zur digitalen Interpolation von Signalen durch Multiplikation verzogerter mit 
einem ersten Taktsignal verkoppelter digitaler Eingangswerte mit entsprechenden Gewichtungsfaktoren (= Filterkoeffi 
z.enten) und Aufadd.erung der multiplizierten Werte. wobei die Gewichtungsfaktoren von einem Zeitdifferenzwert abhan* 
gig sind. der durch den Interpolationszeitpunkt und das Zeitraster des ersten Taktsignals bestimmt wird vergleiche d^n 
Oberbegriff des Anspruchs 1. * l8Mn 

Derartige Verfahren werden bei Interpolationsfiltern verwendet. die insbesondere als Teilschaltungen in digitalen 
Scnaltungssystemen dienen. bei denen eine beliebige Anderung der Abtastrate von digitalen Signalen erforderlich ist 
Systeme. die sich nur mit einfachen ganzzahligen Abtastratenverhaltnissen befassen. sind nichtGegenstand der Erf in ' 
dung. Bei beheb.gen Abtastratenverhaltnissen soil das mit einem ersten Takt abgetastete Signal in eine Abtastfolae 
umgewandelt werden. die mit einem beliebigen zweiten Takt verknupft ist. der im ungOnstigsten Fall sogar ein irrationales 
Frequenzverhaltnis zum ersten Takt aufweist. Auch sollen mit der Erfindung diejenigen Zahlenverhartnisse realisierbar 
sein. bei denen das kfemste gemeinsame Frequenzvielfache beider Takte bei einer unerwonscht hohen Frequenz lieot 
Aus "IEEE. Transactions of Acoustics. Speech and Signal Processing". Band ASSP-32. Nr. 3. Juli 1984 Seiten 577 
^ bis 591 sind unter dem Trtel "Digital Methods for Conversion Between Arbitrary Sampling Frequencies- Autor T. A 
-Ramstad. Verfahren zur beliebigen Anderung von Abtastraten beschrieben. Die zugehdrigen Schaltungen werden als 
hybr.de Systeme bezeichnet. die aus einem ersten Interpolationsfilter mit festem Abtastratenvemaltnis und einem zwei- 
ten Interpolations .Iter (= analoger Resampler) bestehen. Mit dem zweiten Interpolationsfilter werden die Zwischenwerte 
bestimmt. die zertlich beliebig zwischen den festen Abtastwerten des Abtastrasters nach dem zweiten Interpolationsfilter 
liegen und damit beliebige Abtastratenverhaltnisse zulassen. Das erete Interpolationsfilter enthart als Kombination eine 
Interpolationseinrichtung und ein digitales Filter. Mit der Interpolationseinrichtung. die auch als Oberabtasteinrichtuna 
bezeichnet wind, werden -Nuir-Werte entsprechend einem Oberabtastfaktor N zwischen die ursprtinglichen Abtastwerte 
eingefOgt. Mit einer Dezimiereinrichtung werden gegebenenfalls ursprungliche Abtastwerte unterdruckt Erst das nach- 
folgende digitale Filter glattet in beiden Fallen den Verlauf der digitalen Abtastwerte. wobei insbesondere die Sprunqe 
auf die Null-Werte ausgeglichen werden. so daB das Spektrum des Nutzsignals nicht durch Where Frequenzkompo- 
nenten verfalscht wird. Hierzu ist das erste Interpolationsfilter so ausgebildet. daB grOBere FrequenzbereichslQcken in 
dem sich bis ins Unendliche erstreckende Spektrum gebildet werden. Denn auch bei der Oberabtastung gilt daB sich 
die Frequenzspektren bei der halben ursprunglichen Abtastfrequenzen und deren Vielfachen spiegeln Nach'der Inter- 
polationseinrichtung und nach dem digitalen Filter ist allerdings von einer neuen Abtastfrequenz - im Folgenden auch 
erster Abtasttakt genannt - auszugehen. die im ganzzahligen Frequenzverhaltnis zur usprOnglichen Abtastfrequenz 
steht. Das digitale Rlter entfernt dabei wie angegeben die verbleibenden spektralen Kbmponenten zwischen dem Nutz- 
signalband und dem gespiegelten Band bei der neuen Abtastfrequenz und den zugehorigen FrequenzvieHachwerten 
Die Funktionsweise des digitalen Rlters ist dabei einfach diejenige eines digitalen Tiefpasses. der den Nutzsignalbereich 
durchlaBt und die daruberliegenden Frequenzkomponenten unterdruckt. wobei allerdings entsprechend dem Abtast- 
theorem die bekanrrte Spiegelung bei der halben Abtastfrequenz auftritt. Ein digitaler TiefpaB kann also niemals die 
Vielfachen der Abastfreqenz unterdrucken. 

Die spektralen Kbmponenten bei der neuen Abtastfrequenz und den Frequenzvielfechwerten mOssen fur die Rea- 
lisierung beliebiger Abtastratenverhaltnisse durch das zweite. analoge Interpolationsfilter unterdruckt werden. Werden 
diese Storkomponenten nicht unterdruckt. dann werden bei der Erzeugung beliebiger Abtastratenverhaltnisse Storkom- 
ponenten im Nutzsignalband erzeugt. Die Theorie des ersten Interpolations ilters ist hinreichend in "Proceedings of The 
IEEE". Band 61. Nr. 6. Juni 1973. Seiten 692 bis 702 in dem Aufsatz "A Digital Signal Processing Approach to Interpo- 
lation" von R W. Schafer und L R. Rabiner beschrieben. 

Fur das zweite Interpolationsfilter gibt T. A. Ramstad in dem oben genannten Aufsatz in 'IEEE. Transactions On 
Acoustics. Speech and Signal Processing" . Juni 1984. Seiten 577 bis 591 folgende Moglichkerten an: 

1 . Eine Abtasthalteschartung - dieses Rlter hat einfache Nullstellen in der Mitte der noch verbleibenden spektralen 
Komponenten bei den der neuen Abtastfrequenz. Dieser Rltertyp wird auch als einfaches Kammf ilter bezeichnet 

2. Eine lineare Interpolationsschaltung - dieses Filter hatdoppelte Nullstellen in der Mitte der noch verbleibenden 
spektralen Komponenten bei den Frequenzvielfechwerten der neuen Abtastfrequenz. Damit wird eine wirksamere 
UnterdrOckung der Storanteile als bei der einfachen Abtastharteschaltung erreicht. Diese erfordert beispielsweise 
fur em Signal/Rauschverhaitnis von 1 oodB einen Oberabtastfaktor N von 200. der durch das erste Interpolationsfilter 
realisiert werden muB. In der genannten Literaturstelle wild ein -optimales" InterpolationsfiHer angegeben. das eine 
Absenkung des Uberabtastfaktors auf einen Wert von 164 ermoglicht Selbstdie Verwendung sehr rechenaufwen- 
diger Lagrange-lnterpolationsfater hoherer Ordnung erlaubt keine wesentiiche Absenkung der erforderlichen Ober- 
abtastung. Sowohl das lineare Interpolationsfilter als auch das Lagrange-lnterpolationsfilter zeigen ein 
Kammfilterverhalten. 
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Aus EP-A 0 561 067 ist ^fclls ein Verfahren mit einem hybriden SystemRr Abtastratenumsetzung bekannt. 
Dieses System arbeitet allerdings nur mit einem Uberabtastfaktor N = 2 und efreicht damit nur einen relativ schlechtes 
Signal/Rauschverhaitnis. Dies ist jedoch zuiassig. da dieses System fur Videosignalanwendungen gedacht ist. Das 
zweite Interpolationsf ilter ist dabei als TiefpaB realisiert. der al!e Frequenzen oberhalb dem 1 .5-fachen Wert der 

5 ursprunglichen Abtastfrequenz unterdruckt. Das analoge TiefpaBverhalten wird mit einem Transversalfilter erreicht. bei 
dem die Gewichtungsfaktoren der gespeicherten Abtastwerte von einem Zeitdiffererzwert abhangig sind. der im folgen- 
den auch als Inter-Sample-Position bezeichnet wird. Ein derartiger TiefpaG unterdruckt dabei nicht nur die verbleibenden 
spektralen Komponenten bei den Frequenzvielfachwerten der neuen Abtastfrequenz sondern den gesamten Spektral- 
bereich oberhalb seiner Sperrf lanke. Bei einem vergleichbaren DurchlaB/Sperrverhalten ist ein derartiger TiefpaB immer 

10 aufwendiger zu realisieren als eine entsprechende Kammfilteranordnung. Der TiefpaB verarbeitet zwar digitale. abge- 
tastete Werte- in seiner Wirkung stellt er jedoch einen analogen TiefpaB, einen analogen Resampler. dar. 

.Aus "Journal of Audio Engineering Society". Band 41, Nr. 7/8. 1993. Seiten 539 bis 555 ist Von den Autoren R. 
Adams und T Kwan unter dem Tltel "Theory and VLSI Architectures for Asynchronous Sample-Rate Converters" ein 
Verfahren fur ein Abtastratenwandlungssystem bekannt. das einerseits die Verwendung einfacherer Abtasthalteschal- 

is tungen und andererseits die Verwendung von TiefpaBf iltern als analoge Resampler behandelt. FQr einfache Systeme 
werden wegen des geringen Schaltungsaufwandes Abtasthalteschaltungen errpfohlen - fur hoherwertige Systeme wird 
"Njagegen die Verwendung eines TiefpaBf ilters nach dem ersten Interpolationsf ilter vorgeschlagen. das den gesamten 
Frequenzbereich oberhalb des Nutzsignalbereichs unterdruckt. 

Aus den bekannten Uteraturstellen geht hervor. daB nach der N-fachen Uberabtastung und Filterung nach dem 

20 ersten Interpolationsfilter im Frequenzspektrum auf jeden Fall StOrbereiche vorhanden sind. deren Mittenfrequenzen 
bei den Frequenzvielfachwerten der neuen Abtastfrequenz liegen. Die Bandbreite jedes Storbereichs ist dabei gleich 
der doppelten Bandbreite des Nutzsignals. Wenn die Nyquistbedingung fur die ursprungliche Digrtalisierung erfullt ist. 
dann hat die Bandbreite des Storbereichs im Grenzfall maximal den Wert der ursprunglichen Abtastfrequenz. Die Lage 
und Bandbreite samtlicher StOrbereiche ist im Frequenzspektrum durch die ursprungliche Abtastfrequenz und den Uber- 

25 abtastfaktor N def iniert. Die N-fache Uberabtastung der ursprunglichen Abtastfolge bewirkt, daB die relative Breite der 
StOrbereiche bezogen auf die neue Abtastfrequenz urn den Faktor 1/N reduziert wird. Dies erleichtert die Trennung des 
Nutzsignalbandes vom jeweiligen StOrbereich. weil der Ubergangsbereich zwischen DurchlaB- und Sperrbereich fur das 
zweite Interpolationsfilter vergrOBert wird, wodurch der erforderliche Schaltungsaufwand fur das zweite Interpolations- 
filter verringert wird. Dies wird jedoch durch einen hOheren Aufwand fur das Glattungsfilter im ersten Interpolationsfilter 

30 erkauft * 

In der Regel hat man entweder ein sehraufwendiges erstes Interpolationsfilter. z.B. durch eine hone Uberabtastung. 
und ein einfaches zweites Interpolationsfilter. z.B. einen linearen Interpolator, oder aber man hat ein einfaches erstes 
Interpolationsfilter. z.B. mit einer geringen oder gar keinen Uberabtastung. und ein sehr aufwendiges TiefpaBfilter mit 
dem der analoge Resampler realisiert wird. In der letzgenannten Uteraturstelle werden beispielsweise 16.000 feste 

35 Stutzstellen fur das aufwendige TiefpaBfilter gespeichert. 

Aufgabe der Erf indung ist es daher. fur das zweite Interpolationsfilter ein verbessertes Verfahren zur digitalen Inter- 
polation von Signalen anzugeben. welches ein hohes Signal/Rauschverhaitnis ermOglicht. Dies soil mit einem moglichst 
geringen Schaltungsaufwand fur das Gesamtsystem aus erstem und zweitem Interpolationsfilter realisiert werden. - - 
Die LOsung der Aufgabe ergibt sich aus den Merkmalen des Anspruchs 1. Vorteilhafte Weiterbildungen der Erfin- 

40 dungen ergeben sich aus den Unteranspruchen. 

Die Erf indung besteht darin, daB sich die Signaldampfung des zweiten Interpolationsf ilters im wesentlichen kamm- 
f itterartg nur auf die StOrbereiche bei den Frequenzvielfachwerten der neuen Abtastfrequenz beschrankt und fur die 
dazwischenliegenden Bereiche einen beliebigen Signalverlauf zulaBt Dabei weist die Ubertragungsfunktion des zweiten 
Interpolationsf ilters bezuglich jedem der StOrbereiche mindestens zwei Nullstell en auf. wobei die mindestens zwei Nullst- 

45 ellen in jedem StOrbereich nebeneinander liegen oder falls es sich bei den mindestens zwei Nullstellen um eine Nullstelle 
zweiter Ordnung (= zweiten Grades) handelt. mindestens einem der StOrbereiche mindestens eine wertere Nullstelle 
der Ubertragungsfunktion zugeordnet wird. Diese weitere Nullstelle wird vorzugsweise. zur gezielten Verbreiterung des 
Sperrbereiches. neben die vorhandene Nullstelle zweiter Ordnung gesetzt. Durch die mindestens eine weitere Nullstelle 
werden ferner bei den zugehOrigen Frequenzvielfachwerten weitere Nullstellen erzeugt. die dann in der Regel innerhalb 

so der zugehOrigen periodischen StOrbereiche liegen. Dadurch wird die Urrterdruckung der unerwunschten Frequenzbe- 
reich weiter verbessert Zudem nimmt mit steigenden Frequenzen auch in den DurchlaBbereichen des Kammfilters 
dessen Dampfungsverhalten zu. so daB die Lage der Nullstellen innerhalb der StOrbereiche dort immer unkritischer 
wird. Durch die jeweilige Ubertragungsfunktion werden die Gewichtungsfaktoren der miteinander zu verknupfenden 
Abtastwerte im zweiten Interpolationsfilter bestimmt. wobei die Gewichtungssfaktoren wiederum von dem jeweiligen 

55 Zeitdrfferenzwert des Abtastzertpunktes zwischen dem ersten und zweiten Abtastraster abhangig sind. 

Die werteren Nullstellen oder die Modif ikation ihrer Position lassen sich durch folgende Verfahren erreichen. wobei 
die Verfahren miteinander kombiniert>ar und auch mehrfach aufeinanderfolgend anwendbar sind: 

a) die Verwendung einer Abtasthalteschaltung hOherer Ordnung - also mit mindestens dreifachen Nullstellen. 
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b) die Versch.ebung und Addit^^r Impulsantwort im Zeitbereich und Bildung^fc result.erenden neuen Imnuk 
antwort und schlieBlich 'mpuis- 

c) die Verschiebung und Multiplikation der Ubertragungsfunktion im Frequenzbereich und Bildung einer results™ 
den, neuen Impulsantwort. wodurch in jedem Sperrbereich mindesten zwei nebeneinanderliegende Null^iLn 

5 erzeugt werden. y * u "sieiien 

* « 

Die Verfahren nach a) und b) lassen sich mittels einer besonders geeigneten Filterstruktur, der als Name die 
chung "Switchable Time-Continuous FIR" (= STC-FIR) gegeben wird, realisieren. 

Der Zertdifferenzwert bzw. die Irrter-Sample-Position def iniert einen beliebigen Zeitpunkt zwischen zwei festen Abt 
w astzeitpunkten des ersten Frequenzrasters. Durch das erfindungsgemaBe Verfahren laflt sich die weitere bzw lassen 
sich die weiteren Nullstellen bezuglich der Storbereiche beliebig positionieren. so daB eine gezielte StOrunterdrOckuna 
mogl^ch wnd. Durch einen seitlichen Versatz der Nullstellen im jeweiligen Storbereich last sich der jeweilige Damofunas 
bereich einfach in der Breite modifizieren. Durch gezielte Positionierung weiterer Nullstellen kann man beiswelsweise 
den Uberabtastungsfaktor im ersten Interpolationsfilter von N = 200 auf N = 8 reduzieren. wobei das zweitelnteroola. 
'5 t.onsf .Iter nur 5 Abtastwerte miteinander verknCpfen muB. Unter diesen Bedingungen laBt sich ein Signal/Storvemaltnis 
von lOOdB erreichen. Dadurch wird der notwendige Schaftungs- bzw. Programmieraufwand for das Abtastraten- 
-gmwandlungssystem stark reduziert und die Taktfrequenzen bleiben zudem niedrig. 

Die Erfindung und vorteilhafte Ausgestaltungen werden nun anhand der Rguren der Zeichnung naher erlautert: 

?o Fig. 1 zeigt schematisch die grundsatzliche Funktionsweise von Verfahren zur beliebigen Abtastratenumsetzung 
Fig. 2 veranschaulicht die Abtastratenumsetzung an einem Signalabschnitt im Zeitbereich 
Fig. 3 zeigt schematisch die Frequenzspektren. die bei der beliebigen Abtastratenumsetzung auftreten 
Fig. 4 zeigt die Impulsantworten von Abtasthaltefunktionen erster. zweiter. dritter und vierter Ordung gemaB dem 
ersten Verfahren a), 

?s Fig. 5 zeigt die entsprechenden Clbertragungsf unktionen zu den Abtastfunktionen hoherer Ordnung 

Fig. 6 zeigt die Impulsantworten gemaB dem zweiten Verfahren b). das eine Verschiebung im Zeitbereich betrifft 
Fig. 7 zeigt die entsprechenden Clbertragungsfunktionen zu den Zeitverschiebungen ta = T/2 . T/4 und T/8 
Fig. 8 zeigt ein Beispiel for eine resultierende Obertragungsfunktion gemaB dem dritten Verfahren c). das eine 
Verschiebung im Frequenzbereich betrifft, 

»o Fig. 9 zeigt die entsprechende Impulsantwort zu der Frequenzverschiebung, 

Fig. 10 zeigt schematisch als Blockschattbild eine Realisierung des zweiten Interpolationsfilters ahnlich einem FIR- 
Filter, 

Fig. 1 1 zeigt schematisch als Blockschattbild ein vorteilhaftes AusfOhrungsbeispiel des Interpolationsfilters mit einer 
umschaltbaren Matrix. 

» Fig. 1 2 zeigt schematisch im Zeitdiagramm die Oberfuhrung von einem ersten Abtastraster in ein beliebiges zweites 
Abtastraster und 

Fig. 13 zeigt schematisch wie in einer Matrix die verschiedenen Gewichtungsfaktoren in Abhangigkeit von einem 
Zertdiff erenzwert gebikfet werden. 

In Fig. 1 ist als Blockschaltbild ein Schema fur beliebige Abtastratenumsetzungen dargestellt. wie es Oblicherweise 
zur Erlauterung des Verfahrens verwendet wird. Ein analoges Eingangssignal ^ wird mittels eines Analog-Digrtal- 
wandlers 1 digitalisiert, wobei die Digitalisierungsfrequenz f ein ist Selbstverstandlich muB die Nyquist-Bedingung erfullt 
sein. so daB sich der Frequenzbereich des analogen Signals s ain hochstens bis zur halben Abtastfrequenz f oir /2 
erstreckt Nach dem Analog-Digitalwandler kommt eine Interpolationsstufe 2. welche die ureprongliche Abtastrate f 
urn den Faktor N erhoht. Hierzu werden in der Regel "Null"-Werte zu den zusatzlichen Abtastzeitpunkten als Abtastwerte 
eingefGgt Durch die Digitalisienjng mit der Frequenz f ein spiegelt sich das Nutzsignalspektrum bei der halben Abtast- 
frequenz f 8in /2 und wird zwischen den Frequenzwerten f^,, und 1.5-f ein wiedemott und so fort, vgl. Fig. 3. Durch die 
Uberabtastung in einem Interpolator 2.1 mit dem ganzzahligen Faktor N andert sich an dieser periodischen Folge des 
Signalspektrums nichts. Ein digitales Rlter R entfernt jedoch aus den Storspektren moglichst alle Frequenzanteile die 
zwischen der halben ursprunglichen Abtastfrequenz und dem ersten StOrbereich bei der neuen Abtastfrequenz N1 
und zwischen den weiteren SWcbereichen liegen. Damit liegt entsprechend Fig. 3 ein Signal vor. das nur den Frequenz- 
bereich des Nutzsignals enthalt und darOber hinaus nur Storanteile, die bei den VieHachwerten der neuen Abtastfrequenz 
liegen. Wenn dieses Signal mit einem Digrtal-Anatogwandler 3, der mit der Abtastfrequenz fl=N-f ain . arbeitet zuruck- 
gewandelt wird. darm liegt nach einer entsprechenden Glattung angenahert der ursprOngliche Signalverlauf s 9in wieder 
Der Digital/Analog-Wandler3 erzeugt dabei theoretisch nur zu den unendlich kurzen Abtastzeitpunkten der Frequenz 
f 1 Analcgwerte x(t). die durch ein Glattungsfilter h in ein kontinulerliches Signal y(t) umgewandelt weirfen. Das Glat- 
tungsfilter h. das wesentlicher Teil des zweiten Interpolationsfilters 4 ist. kann daher zutreffend auch analoges Interoo- 
lationsf ilter genannt werden. 
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Fur das Verstandnis dSjRtastratenurnwandlung ist davon auszugeher^BiG das Signal y(t) nach dem Giattungs- 
filter h mGglichst dem Eingangssignal s Gin gleich ist. das heiBt daG alle Nutzsignalkomponenten vorhanden sind und 
keine Stflrsignalkomponenten. Ein unterschiedlicher Skalierungsfaktor ist dabei unemeblich. Wenn diese Bedingungen 
ideal erfullt wflren, dann kdnnte mittels einer zweiten Analog-Digital-Wandlereinrichtung 5 das kontinuierliche Signal y(t) 
5 mit einem beliebigen zweiten Abtasttakt f2 digitalisiert werden. Die digrtalisierten Werte wurden erne zweite digrtalisierte 
Abtastfoige y(d) liefern. Wenn das Signal y(d) spater mit einer niederen Abtastfrequenz f aus weiter verarbeitet werden 
soli, dann schlieGt sich der Oigitalisierungsstufe ein digitales TiefpaBfilter Fo und eine Dezimierstufe 6 an, wodurch der 
Abtasttaktf2umdasganzzahlige Frequenzverhaitnis No reduziertwird. Dem Ausgangssignal s aus istdie Abtastfrequenz 
f aus zugeordnet. die sich aus der zweiten Abtastfrequenz f2 mit f ^ = f2/No . ergibt. 
to Aus Fig. 1 ist ohne weiteres erkennbar. daG die Genauigkeit der Abtastratenumwandlung im wesentfichen von der 
Genauigkeit des Signals y(t) abhangig ist. und dieses Signal hangt wiederum von den Eigenschaften des Giattungsf ilters 
h und des digitalen Filters Fi ab. Bei der praktischen Realisierung verzichtet man selbstverstahdlich auf den analogen 
Resampler R in der dargestellten Form, weil er eine zusatzliche Digital/Analog- und Analog/Digital -Umsetzung fur das 
analoge Giattungsf ilter h erfordert. Es sind namlich. wie eingangs beschrieben. auch analoge Giattungsverfahren 
is "bekannt. die mit einzelnen Abtastwerten auskommen und die erforderlichen Zwisch en werte durch Interpolationverfahren 
erzeugen. Den einfachsten Fall stellt hierbei die einfache Abtasthalteschaltung (= latch) dar. die den letzten Abtastwert 
der Folge x(k) so lange festhait. bis der neue Abtastwert vorliegt Eine Verbesserung bringt der lineare Interpolator, der 
die beiden letzten Abtastwerte speichert und daraus durch lineare Interpolation einen Zwischenwert zum gewunschten. 
dazwischenliegenden Interpolationszeitpunkt tP bildet. Andere bekannte Giattungsfilter sind sehr auwendige TiefpaG- 
20 filter Oder das Lagrange-Filter. Da bei alien diesen Rltern die Inter-Sample-Posrtion beliebig wahlbar ist wird digital auf 
diese Weise ein analoges Filter realisiert, daher auch die Bezeichnung "analoger Resampler". 

In Fig. 2 ist vereinfacht das kontinuierliche Signal y(t) uber der Zeitachse t zusammen mit dem ersten und zweiten 
Abtastraster f 1 bzw. f2 dargestellt. Es soil beispielsweise fur den Zeitpunkt tP eine neuer Abtastwert P2 gebildet werden. 
Mit dem analogen Resampler R nach Fig. 1 erfolgt dies einfach durch den zweiten Anaiog/Digital-Wandler 5. der vom 
25 zweiten Abtasttakt f2 zum Zeitpunkt tP getaktet ist. Ohne den Umweg uber den analogen Resampler R muG der Abt- 
astwert P2 jedoch aus benachbarten Abtastwerten P11, P12. P13, P14 der ersten Abtastfoige interpdiert und damit 
berechnetwerden. Derzeitlich unmittelbar vorausgehende Abtastwert P12 undder Interpolationszeitpunkt tP bestimmen 
dabei einen Zeitdifferenzwert At (in der folgenden Beschreibung wird hierfur meist die Bezeichnung "dt" verwendet ) . 
der zur Berechnung des interpotierten Abtastwertes P2 dient. Bei der linearen Interpolation bestimmt er namlich den 
30 Anteilsfaktor aus den beiden benachbarten Abtastwerten P1 2 und P13. Fur die linearen Interpolation sind die weiteren 
Abtastwerte P11 undP14und alle weiteren Werte ohne Interesse. Der Zeitdifferenzwert dtbezieht sich selbstverstand- 
lich auf die Periodendauer T1 des ersten Abtasttaktes f1. In der Zeitdarstellung stellt sich die lineare Interpolation als 
symmetrische, dreiecksfdrmige Impulsantwort h2(t) im Zeitbereich dar. wobei die Symmetrieachse durch den Interpo- 
lationszeitpunkt tP definiert ist. Die zeittiche Breite der Impulsantwort h2(t) umfaGt dabei zwei Abtastperioden T1 und 
35 damit zwei gespeicherte Abtastwerte. Der Anteil- Oder Gewichtungsfaktor ergibt sich aus dem Schnittpunkt der Abtast- 
werte P12 und P13 mit der Impulsantwort h2(t). In Fig. 2 failt der Abtastpunkt P14 es ersten Abtasttaktes f1 mit einem 
Abtastpunkt des zweiten Abtasttaktes f2 zusammen, so daG dort der Zeitdifferenzwert dt zu Null wird. Die Berechnung 
des jeweiligen Anteilfaktors ist in alien Fallen unabhangig von der Frequenz des zweiten Abtasttaktes f2. Lediglich der 
Interpolationszeitpunkt tP ist fur die Berechnung maGgebend. 
<o \n Fig. 3 sind schematisch einige zu Fig. 1 gehOrende Frequenzspektren dargestellt. Das Nutzsignal hat entspre- 
chende dem Nyquist-Kriterium hOchstens die Bandbreite B=f oln / 2 . Dem entspricht eine doppelt so groGe Bandbreite 

Bs bei den StOrbereichen 1 S pS. die mittig zu den FrequenzvieHachwerten des ersten Abtasttaktes f 1 liegen. In der 

ersten Zeile sind gestrichelt die Frequenzspektren angedeutet. die zwischen dem Nutzsignalband und den jeweiligen 
StOrbereichen liegen und die infolge der Spiegelung an der halben ursprunglichen Abtastfrequenz 1 Qin f2 entstehen. 
45 Durch die Oberabtastung mit dem Faktor N andert sich dabei nichts. Die gestrichelten Frequenzbereiche sollen jedoch 
mOglichst gut durch das digitale Filter Fi unterdruckt werden. Das digitale Rlter R ist damit ein typischer digitaler TielpaG. 
dessen Frequenz sich bei der halben Abtastfrequenz f1/2 spiegelt. 

In der zweiten Zeile stellt Fig. 3 den Frequenzverlauf und damit die Ubertragungsfunktion H(f) im Frequenzbereich 
eines mdglichst idealen - und damit analogen - Tielpasses dar. der als Interpolationsfilter 4 alle Storsignale oberhalb 
so des Nutzsignalbandes B unterdruckt. In der darunteiiiegenden Zeile wird der gewOnschte Frequenzverlauf H'(f) eines 
GiattungsfiKers h dargestellt. bei dem kammf ilterartig nach der Erfindung im wesentlichen nur die periodischen StOrbe- 

reiche 1S, 2S pS unterdruckt werden. Der Dampfungsverlauf in den dazwischenliegenden Bereichen ist uner- 

heblich. wenn sie vom digitalen Filter Fi vorher genugend unterdruckt worden sind. In der letzten Zeile wird schematisch 
die spektrale Verteilung y(0 des interpolierten Signals y(t) dargestellt. die bei den einzelnen StSrbereichen nur noch 
55 kleine Signalreste enthait die kleiner als das gewOnschte Signal/Rauschverhaitnis sein sollen. Anstatt der Bezeichung 
H(f) wird insbesondere in den nachfolgenden Gleichungen die gleichwertige Bezeichnung H(a>) verwendet mit 
ca = 2 • it • f . 

Bisher ging man bei derartigen Interpolationsschaltungen von der Vorstellung aus. daG fur hochwertige Systeme 
mOglichst ein hochwertiger TiefpaG als analoger Resampler R verwendet werden muG. der alle Signalanteile oberhalb 
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des Nutzsignalbandes B unterdn^fcsjach dpr PrfinH, J„ ^ ^ 



^ dabei bedeuten: 
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laZu* " Ste, I en * d6St ° 9r6Rer «*■ di * Wmpfun^ Si r ^ U ^' en B da,9este, ' t - Je 9^Ser die Anzah. 
Signal/Rauschverhaitnisses. z.B. bei 80 dB ist eVkennbar SlnE? J * der Betract «ung eines vorgegebenen 
Sperrbereiche verbreitern. Unabhangig 25 JSS^tSXTlS^ Ans " h ' K an *-^d£2JE 
den Frequenzen zu. Am kritischsten "d££»2£^ ^ errt5ereiche in ih ™ Breite m!t steigen 

ledoch nach dem vorbeschriebenen Venahren o2SSJ555?S^ "'^ Fr6qU6nZ6a D '' 6Se kflnnen 

Nach dem erf indungsgemaBen Verfahren tat JUT. z " sa * l,< f ,en Nullstellen versehen werden. 
eine weif ere Nullstelle d^X Ze^e^eSrS ESS?" ,mpu,santwort h W *> *u verSndern. daB mindestens 
Fouriertransformation ergibt sich 962,6,1 P 05 *™ werden kann. Aus dem Verschiebur^aTz '2 



ftt-taJo^e^-Ffo,) 
und durch Addition mil der ursprOnglichen Funktion f(t) folgt: 

,(t) " te > °=> a***" • cos(o>ta/2) • F( w ) 



55 



h(t)-h neo (t) = h a(t (t) + h all (t-ta) 
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Nach diesem VerfahWIrd die resultierende Impulsantwort h(t) = von Fig. 6 aus einer Gberlaaeruna von 

zwe, ze.tverscrK>benen Imputeantworten h2 l( „ und h*(Q eines .inearen .nterpolators h?K-2) JS^SSJS 
Ube^gerungsverfahren.s, eine Addition der beiden ursprunglichen Impu.santworten zur neuen In^u.santTrSi 
Das Abtastfenster der neuen Funkt.cn ist grolier geworden und kann bis zu drei Abtastwerte der eSen^tastfolae^ k) 
s erfassen. Selbstversttndlich k6nnen beliebige Impulsantworten h(t) der gleichen Operation IrJSTSS urn 
ge*eltwe rt ereNustel^ 

^»pSta£^r 2e,tVerSCh0benen ,m P ul ^rt kombiniert werden. um eine zweite zusatz.Se NuSS 

Zum Beispiel sind drei Obertragungsfunktionen einer Abtasthalteschaltung fur verschiedene Zeitverschiebunapn in 
W w 9 - L d l fleSteW - ^"^ientdie eintacne Abtasthalteschaftung ersterOrdung H SH ^S Sner ^StenS^en 
Ver«h,ebung urn ta = tat ergibt sich der mil H SH+ta1 dargesteWe Kurvenverlauf. Natf ^eTzv^te Tzet^en Ver" 
sch,ebung urn ta = ta2 ergibt sich der mit H SH+ta1+ta2 dargestellte Kurvenverlauf. 

Beisp.elswe.se wird mit einer ersten zeitlichen Verschiebung von ta = T1 / 2 die Anzahl aller Nullstellen der unoe- 
c T»1 n t6r ? e " en VerdOPPelt Di6Se StorSte " en entha,ten i e «ie*eils zwei Nullstellen anstatt deTursprOrwTcSn einen 
f 6 m «" r!l n u° Ch ' na,i9en "v**** der zeiUich «" Verschiebung. aber nun mit ta = ?1 / fSewSs eine 
we ere NuHstel.e be, b^immten periodischen StOrstel.en. die bei den geradzah.igen Frequenzvieifachwer eTtf l 6f 1 

S N«e auf " ' d,6Ser SWrbereiche weist damit ™ NullsteHen anstatt der ursprung.ic^n 

20 e .n,2 a p 9 r annte !2 rfahfen C) 2Ur Unterdr0ckun 9 der Storbereiche ergibt sich. wenn in jedem Storbereich 

20 einzelne NuHstelen nebenemnder. .nsbesondere in ftquidistanter Anordnung, durch eine Verschiebung im Freauenz- 
bere.ch posrt.on.ert werden. Dadurch wird die relative Brerte der Oampfungszonen bei den Frequ^zvSfaihweSn von 
M vergrOflert. wodurch die erforderliche Uberabtastung reduziert werden kann. Fig. 13 zeigSzu !m F^eoXdTa" 
gramm e.n Beispiel. Die Ubertragungsfunktion H1 (<o) einer Abtasthaltescharcung erster Ordnung Z d dabe^o rSS 
z.ert dafi e,ne neue Ubertragungsfunktion H neu( o» als Prcdukt von zwe? im Frequenzber^ chTer^hob^en 
25 ursprunglichen Ubertragungsfunktionen H alt (o>) gebildet wird: verscnooenen 

HneuC^) = H2(w) - HI rf ( a - A« ) - HI „( « + A» ) 

ao ST*" rCh c Wird d3S Da i 'T P, " n9SVerha,ten in sdm t"'c"en StOrbereichen so geandert. dafi die jeweiligen Dampfungszonen 
be. den Frequenzv.e.fachwerten von ft relativ breit sind. Obwohl das Durch.aSverha.ten von H2(f) in Fig 8 bTSen 
Frequenzv.etfachwerten von f1 ein relatives Maximum aufweist wirddie StdrsignalunterdrQckung im ersten Storbereich 
m schlechter als etwa - 40 dB. Im zweiten Storbereich liegt die Storsignalumerdruckung bererts unter S dB So I iie 
Stors^nalunterdruokung werter verbessert werden. dann ist es erforderlich die Ubertragungsfunktion H(f) = H3(0 so 
zu modrf.z.erea da6 jeweils eine dritte Nullstelle genau in der jeweiligen Storbereichsmitte positioniert wird: 

35 

H3(to) = H1f» • H2(o>) 

Gegebenenfalls werden eine vierte und funfte Nullstelle zwischen die vorhandenen Nullstellen gesetzt und so weiter 
Dies st mrt dem nachfolgend beschriebenen Venahren moglich. Die resultierende Impulsantwort h2(t) wird durch eine 
Fourier-Rucktransformationwiefolgtgebildet: 

h2(t) = sin(tcol) / o>1 • a(t) + {sin[(2Tl-t) ■ dj - sin(ta>1)} / col • a(t-Tl) + sin[(2T1-t) ■ co1J / a>1 • a(t-2T1) 

piezugehOrigelmpulsantvortrSWzeigtFig. 9. Sie hateinegroSeAhnlichkeit mit der Impulsantwort h(K=2) des linearen 
Interpolators - s.e hat etwa die Form eines gleichschenWigen Dreiecks. aber die beiden Schenkel sind. anders als beim 

J**?!!??'* 0 (K=2) ' ' eiCht konvex 9 ekriimmt - Dj es ist durch die oben beschriebene sinus-Funktion der Impuls- 
antwort bedingt. ^ 

In Fig. 1 0 ist schematise*! als BlockschaltbiW eine Anordnung zur Realisier ung des Interpolate liters 4 dargestellt 
Die Anordnung entspricht bekannten Transversalfilteranordnungen mit zeitabhangigen Koeffizienten. Sie enthait eine 
^eichereir^drtung mr. die eine bestimmte Anzahl 1 bis L bzw. kn bis kn-L + i Abtastwerte des ersten Abtastrasters x{k) 
speichert und d.ese Werte ntittels einer Gewichtungseinrichtung g mit Gewichtungsfaktoren C 0 bis C L . multipliziert und 
in einer Addiereinnchtung ad aufaddiert. 

GleiSung' dUn9 d6r jewel,i9en G^chtungsfaktoren ergibt sich aus der zeitverschobenen Impulsantwort nach folgender 

C ,(dt) = h(l-T1 + dt) mit I = 0. 1, 2 L-1 

Si nt - ei 2* l ? d Z H. den bekannten Vertahren ^rden nach der Erf indung die zeitabhangigen Gewichtungsfaktoren (= 
Koeffizienten) nicht einer gespecherten Tabelle entnommen oder mKtels e^ier umfangreichen Berechnung ermrttelt. 
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abschnrttsweise polynominelle Funktionen liefem ownaoen vei^pn a) und b) mag| IC h. weii diese 

aus einer Speichereinrichtung ausgelesen und mit dem iewpiiin^n ° -1 Fig. 13. werden 

des Interptfationsfilters 4 erweist s?ch dann S£STI1!S?S2^? ^"l*" St ^ ^^beispiel 
nachdem zweiten Veriahren b erzielenSl 2m I SeT^ser! ^ZTTT^ H( ° *" iusat2 «<*e Nullstelle 
P°kto*abschnitts^ 

Gewicmungsfaktoren. In Fig. 1 1 ert4?tfeT£SS" ™ ^^t^ ,0rdie Berech ™«9 der jeweiligen 

Abscbnittsbereiche dt < ta und dt ^^^^^6^^ ^ ^ Matri2en A ' und A " de ™ * 
die abschnittsweise aus geraden StO^Tn1uSnTm e ^e« S^SSTS !T 'T ,santwort "NO dargeste.it. 
NU, AJtastwertevon 1 bis U1. ,hnen **Z£S!Z£££ ^SSS^SS^S * " UmfaSSen 



Abschnitt von 0 < t < ta : 


h(t) = t 


1. 


Abschnitt von ta <t < 1 : 


h(t) = - ta + 2t 


2. 


Abschnitt von 1 < t < 1+ta : 


h(t) = 2-ta 


3. 


Abschnitt von 1+ta < t < 2 : 


h(t)= 4 + ta - 2t 


4. 


Abschnitt von 2 < t < 2+ta : 


h(t) = 2 + ta - 1 


5. 


FOr die nachfoigende Berechnung der Gewichtungsfaktoren geften folgende Gleichungen: 





C ,(di) = h(l + dt) mit dem laufenden Parameter I = 0, 1,2, ... und 
h(t) = h alt (t) + h a(t (t-ta) 

ZwischenstQcke bis zum^ewei igen Sert Z IJSSST'T Beschreibung der geraden 

Impulsantwortenhft). Sie E^ VT^l* ^ di ° Wnf einfachen 

Inpubantwort: 9e ™ e QewcWun ^ aW * en f ^ 

0<dt <* - b Z w. ta<dt<l 



C,(dt) = dt 
C,(dt) = 2 - ta 
C,(dt)= ta-dt 



C,(dt)= -ta + 2dt 
C,(dt)= 2 + ta-2dt 
C,(dt)= 0 
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Dies ergibt fur die Matizer^phd A" folgende umschaltbare KoeffizientenWr Matrixwerte. die als Faktoren fur die 
nachfolgende Multiplikation bzw. als zeitunabhangiger Wert einer Speichereinrichtung entnommen werden: 



5 




0 


1 






A' = 


2-ta 


0 


A" 






ta 


-1 
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-ta 2 
2+ta -2 
0 0 



Fur das dntte Verfahren c) mit der Verschiebung im Frequenzbereich ist die Berechnung der Gewichtungsfakoren nicht 
mehr mrttels einer Matrixanordung wie in Fig. 10 Oder Fig. 1 1 mOglich. sondern sie erfordert eine Berechung von sinus- 
Werten Oder erne Entnahme dieser Werte aus einer gespeicherten sinus-Tabelle. Denn die erforderlichen sinus-Werte 
lassen sich nicht durch eine einfache Polynomreihe approximieren - in diesem Falle wurde man direkt die Polynomreihe 
und damit die MatrixlOsung verwenden. uiynomreine 

^ ,n , 1 2 j st etwas ausfuhrlicher als in Fig. 2 dargestellt. wie aus einer ersten Abtastfolge x(k) eine zweite Abtastfolqe 
y(d) mit beliebigem Frequenzverhaitnis erzeugt wird. Der gesuchte Abtastwert Pd. der zu dem Interpolations- Oder Abt- 
astzeitpunkt tP gehort wird durch eine, insbesondere polynominelle, Interpolation aus der ersten Abtastfolge x(k) qebil- 
det Hierzu wird der Zeitdifferenzwert dt aus dem nachstliegenden Abtastwert Pn der ersten Abtastfolge x(k) und dem 
gewunschten Abtastzeitpunkt tP bestimmt. Entsprechend diesem Zeitdifferenzwert tP werden die Gewichtungsfaktoren 
Ci aus der Ubertragungsfunktion h(t) ermittelt Durch die Speichereinrichtung mr wird dabei sichergestellt da3 die erfor- 
derlichen Abtastwerte Pn und Pn+1 fur die Interpolationsberechnurig von Pd zur Verfugung stehen Die Berechnuna 
der vorausgeheriden und nachfolgenden Abtastwerte Pd-2 bis Pd+2 erfolgt auf ahnliche Weise durch Verschiebunq des 
Interpolationsfensters und Bestimmung des jewdligen Zeitdiffererzwertes dt. 

Patentanspruche 

1 . Verfahren zur digitalen Interpolation von Signalen durch Multiplikation verzOgerter, mit einem ersten Taktsignal (f 1 ) 
verkoppelter Emgangswerte (x(k)) mit entsprechenden Gewichtungsfaktoren (Ci; Q und Aufaddierung der multipli- 
zierten Werte, die von einem Zeitdifferenzwert (At; dt) abhangig sind. der durch den Interpolationszeitpunkt (tP) und 
das Zeitraster des ersten Taktsignals (f1) bestimmt wird, 
qekennzeic hnet durch folgende Merkmate: 

- mittels einer Impulsantwort h(t) im Zeitbereich werden die Gewichtungsfaktoren (Ci; C,) bestimmt wobei die * 
Ubertragungsfunktion H(f ) im Frequenzbereich ein Signaldampfungsverhalten aufweist, das bezuglich der 
Sperrbereiche im wesentlichen auf die bei den Frequenzvieffachen des ersten Taktsignals (f1) liegenden StOr- 
bereiche (1 S. 2S, ... pS, ...) beschrSnkt ist, 

- jedem dieser StOrbereiche werden mindestens zwei Nullstellen der Ubertragungsfunktion H(f) im Frequenzbe- 
reich zugeordnet, wobei 

- die mindestens zwei Nullstellen in jedem Storbereich nebeneinander liegen Oder 

- mindestens einem der Storbereiche (pS) mindestens eine weitere Nullstelle der Ubertragungsfunktion H(f) 
zugeordnet wird, wenn die zwei Nullstellen als Nullstelle zweiter Ordnung (K=2) aufeinander liegen. 

Verfahren nach Anspruch 1 . dadurch gekennzeichnet, daB eine neue Impulsantwort h(t). deren erste weitere Nullst- 
elle der Ubertragungsfunktion H(f) im Frequenzbereich bei einer vorgegebenen Frequenz fn liegen soil, nach fol- 
gender Bemessungsregel aus einer ursprunglichen Impulsantwort h^tO) 

h W = h oeuW = h aftW + h ^(t - ta) mit ta = 1 / 2fn gebildet wird. 

Verfahren nach Anspruch 2. dadurch gekennzeichnet. dafi zusatzliche weitere Nullstellen der Ubertragungsfunktion 
H(f) durch eine wiederholte Anwendung des Verfahrens gemaG Anspruch 2 beliebig positioniert werden kOnnen. 

Verfahren nach Anspruch 1, dadurch gekennzeichnet. daft jedem StOrbereich (1S. 2S pS, ...) mindestens eine 

weitere Nullstelle der Ubertragungsfunktion H(f) im Frequenzbereich zugeordnet ist. so daG alle StOrbereiche min- 
destens drei Nullstellen aufweisen. 

Verfahren nach Anspruch 4, dadurch gekennzeichnet. daB die mindestens drei Nullstellen mittig zu den jeweiligen 
StOrbereichen (1 S pS. ..) bei den Frequenzvielfachwerten ( 1f 1 . 2f 1 pf 1 . . ..) des ersten Taktsignals (f 1 ) liegen. 
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6 BEB^^^^^-^^ 



dabei.bedeuten: 

Tnr die Ordung der mittigen Nullstellen der Gbertragungsfunktion HnfO 

cm bzw. a(t kTt) e,nen Bnhertssprung zum Zeitpun* t bzw. zum verzagerten Zeitpunkt , - kTI. 

H TOU (f) nach folgender Forme! Frequenzwerte + ^fund - A f eine neue Obertragungsfunwlon 

H neu(0 = H a| ,(f-Af)- H^f+Af) 

S- Verfahren nach einem der Anspruche 1 bis 7 riaHnmH na u^ 

Hffl *«, «1C J^sscat^ssrr - " * 

" 10 ' ^^^ZT^^^^^^<^^^^ ( c, c, 

wert (dt) definiert ist. 1 ' °° ^ ^ abschi'ttsweise bezuglich dem Zeitdifferenz- 

Matnxstufe (mx; mxx) gebildet werden. reaiisiert wird. wobe. d.e Gevwchtungsfaktoren mittels einer 
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